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Abstract 

An explicit analytic formula for the transverse and longitudinal im- 
pact factors St,l(N, 7) of the photon using k-r factorization with exact 
gluon kinematics is given. Applications to the QCD dipole model and 
the extraction of the unintegrated gluon structure function from data 
are proposed. 



1 Introduction 

In the present knowledge on perturbative QCD resummations at leading level 
in logarithms of the energy (i.e. logl/xsj) the coupling of external sources, in 
particular a virtual gluon in deep-inelastic reactions, is based on the theorem 
of Jvt factorization |IJ, proven in the leading logarithmic approximation of 
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QCD [|T|. The theorem states that the "unintegrated" gluon distribution, 
i.e. the distribution of energy and transverse momentum of gluons in the 
target, factorizes from the rest of the process. The remaining factor is the 
so-called "impact factor" . Consequently, this "impact factor" is an universal 
quantity, the same in all processes initiated by the same external source, e.g. 
the photon. The "unintegrated" gluon distribution will depend on the target, 
but again the target impact factor can also be factorized out, leaving place to 
an universal interaction term, given by the Balitsky, Fadin, Kuraev Lipatov 
BFKL Pomeron |2|. At next-to-leading level, the modified interaction term 
is now known H] but not yet are the impact factors determined at this order 
of perturbation [[|. It is expected that the effect of the exact kinematics of 
the exchanged gluons, which is our subject, gives the main contribution to 
these higher order terms. 

The fcr factorized impact factors can be conveniently expressed in terms 
of two Mellin variables, 7, conjugated to transverse momentum squared and 
N, conjugated to energy. Up to now, only the impact factors at N = were 
considered, since, strictly speaking, kx factorization has been proven only at 
infinite energy, implying iV = 0. If, however the validity of kx factorization is 
extended to the case of exact gluon kinematics, it implies the knowledge of the 
combined 7, N dependence of impact factors. To our knowledge, this com- 
bined dependence has been derived only for real photoproduction of heavy 
flavors 0. It is the purpose of the present paper to give an explicit analytic 
expression for the 7, N dependence of virtual photon impact factors. To this 
end, using the kx factorization and exact gluon kinematics, we derive the 
explicit formulae for the total cross section of longitudinal and transverse 
photons on any target with a given distribution of gluons. 

In the next section 2, the formulae for the (virtual) photon impact fac- 
tors following from k? factorization are given. The details of the calculation, 
implying multidimensional integration and resummation of generalized hy- 
pergeometric functions are presented in section 3 for the longitudinal case 
and in section 4 for the transverse one. Applications of our results are given 
in section 5.1 for the QCD dipole model || which turned out to be rather suc- 
cessful in description of the deep inelastic total and diffractive cross-section 
of the virtual photons J!]]. In subsection 5.2 we suggest a model independent 
method of extraction from data of the unintegrated gluon structure function. 
Our conclusions are given in the last section. 
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2 An explicit formula for the total cross - 
sections using kr factorization 

Our starting point is the formula for the total longitudinal and transverse 
photon cross-sections given in ||: 

Air 2 ot 

a L = ^rF L = Aaa s Q 2 / dz[z(l - z)} 2 
Q z Jo 

dk 2 r .o / 1 



Using 



we write 



with 



/ k 2 \ 7 ~ 2 

5 L (iV,7) = 4aa s / - 2)]" I dk 2 ' 



2 



l V (Q 2 " A 



d 2 p[ s =, ^ ^ s s (4) 

\p 2 + q 2 {p-ky + Qy {{p-{i-z)ky + Q 2 + k 2 ] N 

where we have used the relation 

(p-(l-z)k) 2 + Q 2 + k 2 } A2 n ._ 2 _ 
x g = x Bj - — ^ Q 2 = z(l-z)Q 2 (5) 

Q 

Similarly, 

cr T = —^-F T = aa s / ofz z + (1 - z) \ 
Q z Jo 

dk 2 r 2 / p p — 



and using (§) and @ we have 



with 



S T (N, 7 ) = aa s £ dz[z* + (1 - zf] J ^ ' 



Q 2 

cPpl^C fcM 2 , , ■ (8) 

It turns out that the integrals (§) and (H) can be explicitly performed and 
expressed in terms of i[> digamma functions. The details of the calculation 
are given in sections 3 and 4. Here we only quote the final formulae: 

q ( i\j \ o vrr(7 + ^+l)r( 7 +l) 
Sl(N,7) = 8aa s 

1 r>( T+ 5)-y,( 7 ) x 3iV 2 -(5 2 -l) J (g) 



(5 2 -l)(5 2 -4) 15 2iV 



ST(J V, 7 ) = 2aa,-^ +i ) r W 



T(N) (6 2 -l)(5 2 -A) 
[ ijj{l+8) - y>( 7 ) iV 2 (3(iV+l) 2 +9)-2iV(^-l) + (g-l)(g-9) 
\ 5 X 4iV 



-i(3(iV+l) 2 +9 + (5 2 -l))} , (10) 



where we adopted the convenient notation 8 = N — 2 7 + 1, 

Note that the poles at 5 = 0, ±1, ±2 in formulae (|9|JT0|) are actually 
absent due to zeroes in the numerators, as it should from regularity of the 
generalized (Meijer) hypergeometric functions appearing in the derivation, 
see later. This provides numerous and non trivial checks of the resummations 
leading to (gMj. 



3 Longitudinal photon impact factor 

As the first step in the calculation we observe that, using the symmetry of 
the integrand with respect to interchange of z and 1 — z, the integrals in (f|) 
can be written as a sum of two terms: 

S L (N, 7) = Saa s (Q 2 ) 2 ~ 7 (A - B) (11) 
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where 



A = l l *w->)rf£ (* 



2 / aK (l„2V 



fp—L, . . (12) 

(p 2 + Q 2 ) 2 [(p- (1 - ^)A;) 2 + Q 2 + 



fl = J £V<i-*)r7£(*T 



^ P ^ , <?>" . . . (13) 

[p 2 + Q 2 }[{p - k) 2 + Q 2 } [(p — (1 - z)k) 2 + Q 2 + k 2 ] N 

Using several times the identity 

1 1 roo 



dt t M ~ l e~ tC (14) 



C M T(M) Jo 
we transform fll~2|) and ([13]) into 



A = TW) H d *[z{l-z)?{&) N J fr {k 2 Yj vdv J dt t N -> 

J d 2 p ex P {-v(p 2 + Q 2 )-t[(p-(l-z)k) 2 + Q 2 + k 2 ]} , (15) 

B =f5firw-')iw/¥ 

dv ( dv' I dt t N ~ l [ d 2 p e -( P 2 44 2 W((^) 2 4<3 2 H[(^(i-)A?) 2 4<3^ 2 ] (16) 



Thus the integration over d 2 p reduces to a gaussian form and can be easily 
performed. After rescaling the variables v = ty, v' = ty' and substitution 
u = tk 2 the integrals over du and dt factorize from the rest and can be 
explicitly evaluated. The final result of these operations reads 

AB irr (7 -l)Tg- 7 + 2) ( ^^ to , (1 _ a) 

hl (y + zr\ (17) 



(1 + 2/) 



JV-27+4 
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dydy'Ky + zKy' + l-z)} 1 -^ 



(l+y + y') N -^ +i 

To evaluate integrals in QT7D one takes h = 1+y, y+z = h[l—(l—z)/h]. Using 
the formula for the Gauss series 

and picking the factor T(7 — 1) in formula (0), the integral over dv leads to 

^ r( 7 -l + ra) (l-z) n 

~ ^ (A^-7 + 2+n)(iV-7 + l+n) n\ 
= 2 Gi(r/-l, JV-7+I; N-j+3; 1-z) , (20) 

where the Meijer function 2 Gi(m, f ; w; t) = Y{u)Y{v)/Y{w) 2Fi(u,v;w;t). 
Substituting this into ([H]) we can integrate over z to obtain 

a - - r(Af -^ + 1) (QT- r„ (2i) 

where 

T A = 3 G 2 (2, 7 -l,iV-7 + l;7 + 3,iV-7 + 3;l) . (22) 



To calculate the integrals over dy and dy' in (|I8 ) we first rescale the variables 
y, y' — > (1 — Then we introduce h = l + y,h' = l + y' and, finally £ = 
ft,//t'. The last change implies 1 + ?/ + ?/ — > 1 + 2?/+ (1 — = — z(l — £)]. 
Using again (|19|) the integral over dydy' can be transformed into the series 



(l+y+y') N - 2 ^ 



1-7 



r*(i - z)]^ ^ N+2 ~^ y r(jv-2 7 +4+n) 

L 1 n NY{N- 2 7 +4) ^ T(iV-7 + 3+n) 12 +U Zj j 

S W - Z) ^ NY% + X% [^(1,^-27+4; iV-7+3;,) 

+ (z->l-z)] . (23) 
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Using known relations [[J on 2 F% functions, one writes also 

dydy'{(y + z)(yi + l-z)] 1 -"< 



(l + y + y') N ~ W 

= ^ 2 " 7 Arr(7 1 _ 1) { 2 G 1 ( 7 -l,iV-7+2;iV- 7 +3^) + (^l-z)} . (24) 
Introducing this into (|T8| ) we perform integration over z to obtain finally 



27rr(iV- 7 + 2)r( 7 + l) 

5 = r^Tv Tb ^ 

where 

T B = (Q 2 y- 2 3 G 2 (3, 7 - 1, N - 7 + 2; 7 + 4, AT - T + 3; 1) . (26) 

It turns out that the series in ( |2"2"| ) and(^) can be summed up and ex- 
pressed in terms of the digamma functions. One uses known relations on 
Meijer functions at z — 1, see ref.|KJ. One writes 

T A = 3 G 2 (2,a,6;a+4,6+2;l) = 
3 G 2 (2,a,6;a+4,6+l; 1) - 3 G 2 (2, a, 6+1; a+4, 6+2; 1) = 

-[(6-1) 3 G 2 (l,a,6-l;a+3 ) 6+l;l)-6 3 G ! 2 (l,a,6;a+3,6+2;l)] . (27) 
and 

T B = 3 G 2 (3, a, 6+1; a + 5, 6+2; 1) = b ^ b ~^ 3 G 2 (1, a, 6- 1; a + 3, 6+2; 1) , (28) 

o 

using the abbreviations 

a = 7-1; 6 = A^+l- 7 . (29) 

Using the generic formula for hypergeometric functions 3 -F 2 (l, a,, 6; a+p+ 
1, 6+5+I; 1) iTT|, one writes: 

3 G 2 (l,a,6;a+p+l,6+g+l; 1) = 
( _ 1)pH T(p+q+l) T(b-a-p) _ 

1 ij r( P +i)r(g+i)r(6-a+g+i) ma) nb)) + 

T(b-a-p) ^ 1 1 T(l + q+k)T(l-a-p+k) 
} r(l-a)r(g+l) ^ k-pT{k+l) T(b-a+q-p+k + l) 

+ {a <-> b,p <-> q}} . (30) 
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All in all, using for convenience the notations 5 = b — a — 1 and N = b + a, 

we get 



2T B _ 1 f ^ a + i + S)-^(g + l) 3N 2 -(5 2 -l) ^ 

JA iV " (52_i )(5 2_ 4) ^ 5 27V 

4 Transverse photon impact factor 

ot can be evaluated along the similar lines as a^. Here we mark only the 
main differences. The key point is that a? can be evaluated using cr^, which 
gives a noticeable simplification of the painful calculation. 

Let us come back to the calculation of cr^, starting with formula ( pT|) and 
define 

(DL) = Q 2 A ; (NDL) = Q 2 B . (32) 

Using (ffi) we write 



S T {-y,N) = 2aa s ((DT) - (NDT)) (Q 2 f 7 (33) 



with 



(DT) = £ dz[z 2 + (1 - zf] J dk 2 (k 2 y~ 2 



,e P ^t {Q y , . (34) 



(p 2 + Q 2 ) 2 [(p-(l-z)k) 2 + Q 2 + k 

(NDT) = J 1 dz[z 2 + (1 - z) 2 ] J dk 2 (k 2 y- 2 
A ggdj) (35) 

[p 2 + g 2 ][G? - A;) 2 + g 2 ] [(p - (i - z )kf + g 2 + fc2]iv 

Let us rewrite the quantities of interest in the following form: 

(DL) = f 1 dzA L (z) I(z) 
Jo 

(NDL) = ( dzAAz) J(z) 
Jo 

(DT) = ( l dzA T (z) [K(z)-I(z)\ 
Jo 

(NDT) = f 1 dzA T (z) \K(z) - J(z) - \l(z)} , (36) 
Jo 2 



where 

A L {z) = z(l - z) ■ A T (z) = z 2 + (l-z) 2 = l- 2A L {z) . (37) 
The integrals /, J, K, L are defined as follows: 

I(z) = I ^ [ d 2 v~^ (38) 

J k A J F (j ) 2 + Q2)2 (1 _ z )ky + Q 2 + k 2 ] N ' 



J(z)=r d ^ [d% & {Q2)N (39) 

J ^ J F i P 2 + Q2]i(j)-ky+&}[(p-(i-z)ky+Q 2 VkT 

K( Z )=[ d 4^r^ P ^ , , , (40) 

J J F (pl + Q 2 ) [(p- (1 _ 2)fc)2 + Q2 + J.2JJV 



= f—k 2 ^ f d 2 v — (41) 

^ ' J ^ J F \p*+Q*][(p-k) 2 + Q 2 ] ^_^_ Z ^2 + Q2 + ^N ' 

Consequently, one reads for (j L : 

(£>L) - (NDL) = f 1 A L (z)[I(z) - J{z)\ . (42) 
Jo 

The corresponding expression for cry shows a nice simplification, since the 
integrand K(z) cancels from the calculation: 

(DT) - (NDT) = £ A T {z)[J{z) - I{z) + l -L(z)] . (43) 
Using (p7|), we get 

(DT) - (NDT) = 2[(DL) - (NDL)] + 

+ f 1 dz[J(z) - I(z)] + - / dz(l - 2z(l - z))L(z) . (44) 
Jo 2 Jo 



From formula (20) , it is straightforward to obtain (in the a, b notations): 



, 2a r(fo + i) a 



H?) = Q Y(bTtf a 2Gl(a ' b] b + 2; 1 " z) ' (45) 
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where 2 G?i(a, b; b + 2; 1 - z) = r(o)r(6)/r(6 + 2) 2 Fi(a, 6; b + 2; 1 - z). 

After some straightforward transformations from formula (|25|), one writes 

J(^)=Q 2a r ,y + | ) 1 ^ (l-^)[^ 1 2 G 1 (a,6+l;6+2;l-^) + (^ 1-^)] , (46) 
r(0+a + l) 

L ^ = Q 2a m+JT 2Gl{a + h b] b + 1; 1 " z) + {z ^ 1 ~ z)] ' (47) 

where is obtained from J(z) by suppressing the factor z(l — z) and 
replacing a — > a + l,b — > b — 1, due to the numerator fc 2 instead of Q 2 for 
J(z) in formula (41). 

Inserting expressions ([|f|[|6|,f|7]) in formula fl44]) , one gets after integration 
over z: 

(DT) - (NDT) = 2[(DL) - (NDL)] + ^ Q 2 " 

x [ 3 G 2 (1, a + 1, 6; a + 2, 6 + 1; 1) - ba 3 G 2 (l, a, 6; a + 2, b + 2; 1) 

-(6-l)(a + l) 3 G 2 (l,o + l,6-l;a + 3,6 + l;l)] . (48) 

All in all, using for convenience the notations 5 = b — a and N = b + a, one 
finally obtains 

(DT) - (NDT) = r( 5 (a +V 1 



T(a+6) ^ (5 2 -l)(5 2 -4) 
' jj(a+l + 5)-jj(a+l) A^ 2 (3(A^+l) 2 + 9) -2iV(<5 2 -l) + (£ 2 -l)(£ 2 -9) 

_ X 



4iV 

2 



-i(3(iV+l) 2 +9+(5 2 -l))} (49) 



5 Applications 

5.1 Comparison with the dipole model 

In a recent paper [|l2j we have shown that the dipole model is, strictly speak- 
ing, not compatible with the formulae obtained from the k^ factorization 
including the exact kinematics of the corresponding Feynman diagrams. The 
point is that the results from kx factorization are non-diagonal in impact pa- 
rameter space, contrary to the fundamental assumption of the dipole model. 
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Thus, it is worthwhile to derive the explicit modifications of the model due 
to the gluon kinematics. 

The formulae of the previous section can now be compared with those 
obtained in the QCD dipole model, which amounts to consider [^] the impact 
factors at N — 0.: 



Q di P( ,_ Q(AT n v 7T 2 aa s 7(1-7) r 2 (i- 7 )r 2 ( 7 ) 

(7) = ^=o, 7 ) = 3 x _ |7 r(l-7)r(|+7) (50) 



s*( 7 ) = ^=o l7 ) = ^ (i+7)( ^^r 3 (i 7i 2(7) 



7 r(|-7)r(|+ T ) ' 



(51) 



The knowledge of Sl 7 t(N, 7) allows one to take into account the modifi- 
cations of the QCD dipole model due to exact gluon kinematics. Indeed, let 
us insert the BFKL pole in the formula for the unintegrated gluon structure 
function 

9n(i) = ^l^ (go)" 7 ;x(7) = 2^(l)-#y)-lKl-7), (52) 

where a is the effective value of the strong coupling constant a = in 
the BFKL kernel, ^(7) is the Mellin transform of the probability to find a 
dipole in the target, Q sets the typical model scale of the unintegrated gluon 
structure function g(x g ,k 2 ), (see formula (0)) and 

V{1) = 2^WTT) (53) 

is the Mellin transform of the probability to find a gluon in a dipole. In its 
final form, the modified QCD dipole model for longitudinal (||]) and transverse 
(0) cross sections with full kinematics read: 

°l = Q~ 2 j ^- (x Bl r^MiMi)s L (N= ax( 7 ),7) (U) 7 , (54) 

°t = Q- 2 J ^-{x B] r & ^v{ 1 )w{ 1 )S T {N= a X (7),7) (j^Y ■ (55) 
Note that in these formulae, one has 5 = N — 27+1 = axil) — 27 + 1. 



11 



Formulae ( Eg ) and give the explicit dependence of the cross-sections 
in the shift of the hard pomeron intercept. 

We have compared numerical results (in the saddle-point approximation) 
from the Eqs (0) and (|T0|) with those obtained from Eqs (|50| ) and (|5lD in the 



range .01 > x > .0001 and 20Gev 2 < Q 2 < 160Gev 2 . It turns out that they 
give a rather similar dependence on both x and Q 2 (the deviations do not 
exceed 5%). Normalization changes, however: the cross-sections including 
the full gluon kinematics are by about factor 2 smaller than the ones obtained 
from the high-energy approximation. Note the interesting fact that the ratio 
R = <7l/<Jt is affected: it increases by about 15%. 

It is not very surprizing that the dependence on kinematic variables does 
not substantially differ in the two approaches. Indeed, this dependence is 
mostly controlled by the position of the saddle point which is the same in 
the two formulae. The rather important change in normalization, however, 
could not have been easily guessed. In particular the ratio of normalizations 
in R = <jl/<jt is a quantity of experimental interest. 

5.2 Method of extraction of the unintegrated gluon 
distribution 

The knowledge of the impact factors Sl(N,j) and St(N,j) as explicit ana- 
lytic functions of their two variables allows a model independent determina- 
tion of the unintegrated gluon structure function from or ot or from the 
experimental observable F 2 = -£^(&t + 
Indeed, formulae (^JT]) yield 

Zt,l(N,j) = [^L(x Bj f fdQ 2 (Q 2 y\ T!L (x Bj ,Q 2 ) 

J Ob Q j J 

= 9N(i)S T ,L(N,l) ■ (56) 

Hence, from formula (fj) defining the unintegrated gluon structure function, 
one gets 

(57) 



This relation (57) has quite interesting features. 

First, it allows a determination of the unintegrated gluon structure func- 
tion from experiments. For instance one may consider the ratio (<7t + 
pl)/{St + Sl) using data on F 2 . 
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Second, the expected universality properties of g(x g , k 2 ) from kx fac- 
torization gives predictions for various processes e.g. the ratio R = 
the photoproduction and leptoproduction of heavy flavors, diffractive lepto- 
production of vector mesons and any other process where kx factorization 
applies. 

Third, the applicability of this relation goes beyond a specific model like 
the QCD dipole one, provided the coupling to the virtual photon comes from 
the exchange of two off-mass shell gluons with exact kinematics. 

These results comes from the assumption that fey factorization is a valid 
approximation, e.g. the dominance of the two considered Feynman diagrams 
for the impact factors, when exact gluon kinematics can be considered. It 
remains to know to what extent higher order QCD contributions and non 
perturbative corrections may spoil this approximation. 

6 Conclusions 

Our formulae (|9|,|l0D give the exact two- variable dependence of the longitudi- 
nal and transverse impact factors of the photon in terms of the usual Mellin 
variables N, 7. N is the variable conjugated to XBj, while 7 is conjugated to 
Q 2 . In the particular framework of the QCD dipole model, it gives the mod- 
ification taking into account the shift in N = ax{l) of the BFKL singularity. 
This results mostly in a change of the normalization of the cross-section, 
whereas the dependence on kinematic variables predicted by the QCD dipole 
model is hardly affected. Note that the relative normalization in R = ^ is 
increased by 15%. 

More generally, our result opens the way of a model-independent extrac- 
tion of the universal unintegrated gluon structure function which appears in 
various processes, whenever kx factorization can be applied. 
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